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In this paper we revisit non-spherical perturbations of the Schwarzschild black hole in the context 
of f(R) gravity. Previous studies were able to demonstrate the stability of the f(R) Schwarzschild 
black hole against gravitational perturbations in both the even and odd parity sectors. In particular, 
it was seen that the Regge-Wheeler and Zerilli equations in f(R) gravity obey the same equations 
as their General Relativity counterparts. More recently, the 1+1+2 semi-tetrad formalism has been 
used to derive a set of two wave equations: one for transverse, trace-free (tensor) perturbations and 
one for the additional scalar modes that characterise fourth-order theories of gravitation. The master 
variable governing tensor perturbations was shown to be a modified Regge-Wheeler tensor obeying 
the same equation as in General Relativity. However, it is well known that there is a non-uniqueness 
in the definition of the master variable. In this paper we derive a set of two perturbation variables 
and their concomitant wave equations that describe gravitational perturbations in a covariant and 
gauge invariant manner. These variables can be related to the Newman-Penrose (NP) Weyl scalars 
as well as the master variables from the 2+2 formalism. As a byproduct of this study, we also 
derive a set of useful results relating the NP formalism to the 1+1+2 formalism valid for LRS-II 
spacetimes. 


I. Introduction 

General Relativity (GR) has been a hugely successful theory enabling us to explain gravitational phenomena on both 
astrophysical and cosmological scales. The theory has passed many rigorous tests, such as solar system constraints 
[55] or tests using binary pulsars [50]. However, these tests only involve weak gravitational fields and/or velocities 
v c, meaning that General Relativity is still essentially untested in the strong-field, highly dynamical v ~ c 
regime, where high-energy corrections to gravitation may occur. Examples of such strong field systems include black 
holes and neutron stars, which are of particular interest for existing (Advanced LIGO/Virgo) [26] and future GW 
detectors. These systems are also expected to provide a testing ground in which we can understand the dynamical and 
phenomenological features of modified theories of gravitation in the strong field regime. In this paper, we revisit linear 
perturbations to the Schwarzschild black hole in both GR and f(R) theories using the 1+1+2 semi-tetrad formalism. 
There are a few advantages to this approach. The first is that the system of equations describing the structure of the 
spacetime will be a set of first order differential equations in the physical curvature and dynamical variables of the 
covariant derivatives of tetrad vectors, unlike metric approaches which express the Einstein field equations (EFE) as 
a system of second order partial differential equations in terms of metric functions. Secondly, gauge-invariance will 
be naturally realised via the Stewart-Walker lemma [48]. Finally, all objects will have a well defined physical and 
geometrical meaning, as the variables are related to the kinematical properties of timelike and spacelike congruences. 

Some of the most natural extensions to General Relativity are those that appear in the low energy limit of various 
fundamental theories. Among the most popular candidates for ultraviolet modifications to GR are the set of fourth 
order theories of gravity (FOG), for which the Einstein-Hilbert action is modified by additional terms that lead to a 
set of field equations that are fourth order in the metric tensor. A particular subclass of fourth order theories that 
has received a lot of attention in the literature are the f(R) theories in which the modification is some function of the 
Ricci scalar. A motivation for considering such a class of theories can be taken from the equivalence of metric-/(l?) 
gravity to a scalar-tensor (ST) counterpart, in which the gravitational interaction is mediated by the spin-2 graviton 
degrees of freedom as well as a non-minimally coupled spin-0 scalar degree of freedom. From a fundamental point of 
view, such ST theories arise as a natural by-product of string theory (e.g. [11, 12, 14]) and date back to the original 
work by Jordan [27], Fierz [19] and Brans and Dicke [5]. 

In order for f(R) theories to be viable alternatives to General Relativity, there are a number of minimal con¬ 
straints that we can impose. For example, they must reproduce cosmological dynamics consistent with observations, 
they must be free from tachyonic instabilities and ghosts and they must reproduce acceptable Newtonian and post- 
Newtonian limits. In addition, we can also demand that certain well defined solutions in General Relativity, such as 
the Schwarzschild solution, be stable against perturbations in any theory of modified gravity. The analysis of linear 
non-spherical perturbations of a Schwarzschild black hole in ST theories has previously been considered in [15, 35, 41]. 

In General Relativity, the Jebsen-Birkhoff theorem constrains spherically symmetric vacuum spacetimes to be 
either static or spatially homogeneous. It has been shown that the rigidity of the theorem is upheld even upon the 
introduction of perturbations. Notably, almost spherical symmetry and/or almost vacuum imply an almost static or 
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almost spatially homogeneous spacetime [18, 23, 24]. However, it is not clear that this theorem should necessarily hold 
for theories with additional degrees of freedom, as is the case for scalar-tensor and fourth order theories. Importantly, 
[40] found a non-zero measure in the parameter space of f{R) theories for which a Jebsen-Birkhoff like theorem exists. 
This provides us with an additional set of constraints that guarantee the stability of a Schwarzschild solution under 
generic perturbations. Using these results, [41] applied the 1+1+2 semi-tetrad formalism to study non-spherical linear 
perturbations to the Schwarzschild black hole in f{R) gravity. Working in the Jordan frame, [41] derived a modified 
Regge-Wheeler tensor that unifies the axial and polar degrees of freedom into a single transverse-traceless tensor 
obeying a tensorial form of the Regge-Wheeler equation [9, 43]. 

In this paper, however, we advocate that a cleaner way to study the evolution of gravitational perturbations in the 
Schwarzschild spacetime is via the Weyl tensor. In particular, we introduce a perturbation variable J a b constructed 
from the electric and magnetic Weyl 2-tensors. We show that this variable obeys a covariant, gauge-invariant and 
frame-invariant wave equation that is valid in both the polar and axial sectors. As such, this variable unifies the 
Regge-Wheeler and Zerilli equations into a single compact form. Unlike the modified Regge-Wheeler variable for 
f(R) gravity, this perturbation variable demonstrates the decoupling of the scalar and gravitational wave modes in 
a clean and transparent manner. All equations presented here will be valid in both f(R) and General Relativity 
unless explicitly stated otherwise. Secondly, we introduce a purely axial variable V a b that constitutes a covariant, 
gauge-invariant generalisation of the Regge-Wheeler (RW) variable. This variable can also be used to re-interpret the 
RW term as the radial part of the magnetic Weyl tensor. The variable V a b also obeys a covariant and gauge-invariant 
wave equation that reduces to the RW equation when restricting ourselves to the RW tortoise coordinates. 

Finally, we argue that the 1+1+2 formalism provides a powerful tool for understanding the physical and geometrical 
interpretation of other approaches to perturbation theory in General Relativity. Notably, we see that the Newman- 
Penrose (NP) Weyl scalars 4 ' 0 and +4 can be related to our perturbation variable J ab and the imaginary part of the 
NP Weyl scalar 4+ can be related to our perturbation variable V a b ■ Likewise, for the 2+2 formalism, the polar master 
variables x anc l + are just related to the electric Weyl 2-tensor £ ab and the master variable ? is related to R a b■ In 
the axial sector, the master variable n can be exactly related to V a b , as is expected. A roadmap to the set of master 
equations in the 2+2 formalism for f{R) gravity is highlighted based on the results presented here. 

This paper is organised as follows. In Sec. II we summarise fourth order theories of gravity and summarise the 
approach taken to model the f{R) theory. In Sec. Ill we introduce the 1+3 formalism as the starting point for the 
1+1+2 approach used here. Sec. IV builds on the previous Section and introduces the 1+1+2 formalism and defines 
all necessary variables and operators that will be required in subsequent sections. In Sec. V we start with review 
of the previous approaches to studying gravitational perturbations in the Schwarzschild spacetime with the 1+1+2 
approach including the extension to f(R)- We then introduce a set of new perturbation variables derived from the 
Weyl tensor that are central to the remainder of the paper. In Sec. VII we use the results from the previous section 
to provide intuition and insight into other approaches to perturbations of the Schwarzschild spacetime. Notably, we 
discuss the geometric interpretation of the NP Weyl and NP Ricci scalars before discussing the physical interpretation 
of the 2+2 gauge invariant master variables. 


II. Modified Gravity and Black Holes 
A. Fourth Order Gravity 

General Relativity is a unique four-dimensional theory for which the gravitational interactions are mediated by a 
massless spin-2 particle, the graviton, and the field equations are second order. In fourth order gravity, we consider 
the following modification to the Einstein-Hilbert action 

J d A x R -+ J d 4 xV=g f (R,RabR ab ,C abcd C abcd ) . ( 2 . 1 ) 

As the curvature functions contain second derivatives of the metric, the resulting field equations are fourth order. 
In fact, a consequence of Lovelock’s theorem is that the field equations for a metric theory of modified gravity in a 
four-dimensional Riemannian manifold will admit higher than second order derivatives [29, 30]. Such higher order 
derivatives are potentially problematic as they typically lead to Ostrogradski instabilities [42]. The sub-class of f(R) 
theories evade these instabilities as they are degenerate, i.e. the highest derivative term cannot be written as a function 
of canonical variables. This can be seen through the fact that in the Ricci scalar R , only a single component of the 
metric appears with second derivatives. The concomitant new degree of freedom can then be completely fixed by the 
goo constraint preventing ghost instabilities from arising in f(R) theories. Alternatively, the f(R) class of theories 
can simply be re-cast as a scalar-tensor theory for which the gravitational interaction is mediated by a spin -0 scalar 
as well as the spin-2 graviton degrees of freedom [46]. Were we to consider theories involving curvature invariants 
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R a bR ab , R a bcdR abcd or Weyl invariants C ab cdC abcd , these would be non-degenerate and our resulting theory would be 
plagued by Ostrogradski instabilities. In this paper, we restrict ourselves to the Ostrogradski stable f(R) theories. 


B. f(R) Gravity 


The action for f(R) gravity can be written as a simple generalisation of the Einstein-Hilbert action 

s = \ J dV f( R ) + ^M(gab^)] , 


( 2 . 2 ) 


where Cm corresponds to matter fields present in the theory. Upon variation with respect to the metric, the perturbed 
action can be written as 


SS = - 


1 

2 


dV f^g 


\fg a b5g' 


f'SR + T^l Sg ab 


(2.3) 


where, following the notation of [41], a ’ denotes differentiation with respect to R and is the standard energy- 
momentum tensor of matter. Demanding that the action be stationary with respect to variations in the metric, it 
follows that 

f (Rab - \g ab R} = \gab (/ - Rf) + v 0 v 6 /' - g ab nf + t™. ( 2 . 4 ) 


Clearly, if f = R then the field equations reduce to those of General Relativity. In the covariant approach to fourth 
order gravity, the field equations for f(R) can be re-expressed in the form of effective Einstein field equations 


G ab = 




= T, 


ab 


(2.5) 


where T ab is now a combination of two effective fluids corresponding to an effective curvature fluid and an effective 
matter fluid, 


rp _ rpM I rpR 

-L ab J-ab ' ab 

rpM 

rpM _ ^ ab 

as- f , 


rpR _ ±_ 
1 ab ft 


\gab (/ - Rf) + VaVbf 


gabOf 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


The key advantage to this representation is that it is much easier to adapt techniques from the covariant approaches 
to relativistic cosmology to f(R) gravity [17, 41]. 

In the field equations presented above, we have terms involving higher than second order derivatives of the metric. 
For example, in Eq. (2.8), we have terms of the form S7 a S7bf ■ I n non-degenerate theories, such terms would be 
problematic. 

As per [41], the calculations are restricted to the Jordan frame for which matter is minimally coupled and the 
gravitational scalar is non-minimally coupled to the curvature. In this frame, the dynamics of the extra gravitational 
degree of freedom will be determined by trace of the effective Einstein field equations (EFE) resulting in a linearised 
scalar wave equation for the Ricci scalar. 


C. Schwarzschild Black Hole in f(R) 

In [39, 40], it was demonstrated that for f(R) theories, where / is of class C 3 at R = 0 with the conditions /(0) = 0 
and /'(0) ^ 0, any almost spherically symmetric solution with almost vanishing Ricci scalar in empty space for some 
open set S will be locally almost equivalent a part of the maximally extended Schwarzschild solution in S. This 
important result constitutes a Jebsen-Birkhoff like theorem that details the conditions required for the existence of 
Schwarzschild spacetimes in f(R) theories. 

As / is of class C 3 at R = 0, it is implied that [40] 

imi < + 00 , |/"(0)| < +oo, |/"'(0)| < +oo, 


(2.9) 
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with the constraint that [40] 


/(0) = 0 R = 0. 

Performing a Taylor series expansion about the background, the linear order result yields 

m = foR- 


( 2 . 10 ) 


( 2 . 11 ) 


III. The 1+3 Covariant Formalism 
A. The Formalism 

We use a covariant and gauge-invariant approach built on a local 1+3 threading of spacetime with respect to a 
preferred timelike congruence [17, 51]. Consider a spacetime (A i,g) and introduce a preferred timelike congruence u a 
generated by a set of observers with 4-velocity 

dqr a 

U a = u a u a = - 1, (3.1) 

where r is the proper time as measured by the fundamental observers. The spacetime is then locally split in the form 
R® H with R denoting the timeline along u a and H denoting the tangent 3-spaces perpendicular to u a . The induced 
metric tensor on the 3-spaces is h a b = gab + u a Ub and can be used to project tensors onto the spatial surfaces. For 
example, any 4-vector X a can be decomposed into a 3-vector component defined on the 3-space and a (1+3) scalar 
component parallel to the congruence: 

X a = Xf + u a X || where Xf = h a b X b and X\\ = u b X b (3.2) 

Similarly, any second rank tensor may be irreducibly split into scalar, 3-vector and projected, symmetric, trace-free 
(PSTF) 3-tensor components: 


X 


ab 



scalar 


3-vector 




3-tensor 


(3.3) 


where the 3-vector component can be re-expressed in terms of a genuine 3-vector and the alternating tensor e a b c = 
u d rjdabc- Together, the timelike congruence u a and induced metric h ab allow us to define two preferred derivative 
operators. The first operator is a convective derivative along u a 

X a ... b c - d :=u e X e X a ^- d (3.4) 


and the second is a totally projected spatial derivative defined in the 3-surfaces orthogonal to u a 


De -w-':,..,/ := I>"„, •. • h“„h,“ ■ ■ ■ Kh\ V,,W" 


(3.5) 


B. The Variables 

The three main groups of variables in the 1+3 formalism are the set of kinematical variables describing the flow of 
u a , the set of variables associated with the gravitational field as well as the variables describing the energy-momentum 
tensor. The kinematical variables arise from the irreducible decomposition of the covariant derivative of the timelike 
congruence: 


VfcU a = U a U b + -Oh ab + (Tab + UJab, (3.6) 

where u a is the acceleration of the congruence, a ab the volume preserving shear, oj ab the vorticity and 6 the expansion 
of the congruence. These variables describe the kinematical properties of the flow associated to the observer u a . 

The gravitational field is covariantly described by the Ricci tensor and the Weyl tensor. The Ricci tensor describes 
the local gravitational field at a given point due to the local matter content, as can be seen via the EFE. The Weyl 
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tensor describes the non-local gravitational field that is mediated by gravitational waves and tidal forces. These two 
objects allow us to decompose the gravitational field into the local and non-local terms as follows: 

1 R 

Rabcd — C a bcd T ~ ( Qac^bd T 9bdRac QbcRad QadRbc) 7 {^QacQbd QadQbc) • (3.7) 

2 0 

In addition, we can decompose the Weyl tensor into its irreducible parts 

E a b = C acbd u c u d and H ab = ^e a cd C cdbe u e , (3.8) 

where E ab denotes the electric Weyl tensor and H ab the magnetic Weyl tensor. The electric Weyl tensor generalises 
the Newtonian tidal tensor whereas the magnetic Weyl tensor has no Newtonian counterpart as it describes genuine 
relativistic effects such as frame dragging. 


C. The Equations 

In the 1+3 approach, the evolution and constraint equations for the variables arises from two sets of identities 
satisfied by the Riemann tensor with the EFE only being used to algebraically substitute the Ricci tensor for the 
equivalent energy-momentum terms. The first set of equations are given by the Ricci identities applied to u a 


— E abcd U , 


(3.9) 


a second set of equations are provided by the once contracted Bianchi identities, 

^ C a bcd ^[6^a]c T ^9c[b^a]E> 

and the last set of equations are provided by the twice contracted Bianchi identities 

V b Tab = 0. 

The algebraic relations for the Ricci tensor are given by EFE 

Tab Tab Tcjab T A Q a b • 


(3.10) 


(3.11) 


(3.12) 


IV. The 1 + 1+2 Covariant Formalism 
A. The Formalism 

The 1+3 formalism has been hugely successful in relativistic cosmology as it is particularly adapted to the models 
that it aims to describe. For homogeneous and isotropic cosmologies, the 1+3 formalism covariantly factorises out the 
only essential coordinate, time. The background spacetime is described by a system of ODEs for a set of 1+3 covariant 
scalar variables and all conventional ODE methods may be invoked. As the 3-vectors and 3-tensors vanish in the 
background due to the symmetry, they are implicitly gauge invariant as per the Stewart-Walker lemma [48] . Covariant 
and gauge-invariant linear is realised in this approach by noting that the non-vanishing background variables will be 
zeroth order and all variables that vanish in the background will be first order. We can then linearise the full system 
of equations about a specified background, dropping all terms that are second order and higher. Under this procedure, 
all vector-tensor and tensor-tensor couplings are killed off as they are implicitly second order or higher. The result is 
a system of equations that is tractable and may be solved. 

If we relax the symmetry of the spacetime, there will be non-zero vectors and tensors in the background and the 
1+3 system of equations is rendered intractable. However, if the spacetime possesses a preferred spatial direction at 
each point (such as spherical symmetry, local rotational symmetry or G 2 spacetimes) we can introduce an additional 
frame vector such that it covariantly factors out the two essential coordinates: time and the preferred radial direction. 
Consequently, all vectors and tensors will vanish on the 2-surfaces of homogeneity resulting in a system of ODEs 
for 1+1+2 scalar variables. Historically, such decompositions were originally introduced by [25] and were further 
developed by [4, 7-9, 33, 52, 53]. In this paper, we follow the approach to the 1+1+2 formalism as given in [8, 9]. 
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In the 1+1+2 formalism we start by introducing a preferred spacelikc congruence n a that is orthogonal to u a such 
that 


n a u a = 0, n a n a = 1. (4.1) 

The induced metric on the tangent 2-surfaces orthogonal to n a and u a is defined by 

Nab = K b - n a n b = g ab + u a u b - n a n b , (4.2) 

and can be used to project all vectors and tensors orthogonal to u a : n a N ab = u a N ab = 0. Any 3-vector tp a may be 
irreducibly split into a scalar parallel to n a and a 2-vector \& a lying in the 2-surfaces orthogonal to n a : 

ipa = *ri a + Tv (4.3) 

Similarly, any 3-tensor ip ab can be irreducibly split into a scalar, 2-vector and PSTF (with respect to n a ) 2-tensor: 


where 


1 


2 -vector 


i’ab = (j'ianfc “ ~^ Nah j + 2 ^(a™6) + ^{afc} , 

scalar ^ 

2- tensor 


* = n a n b ip ab , * a = N a b n c ip bc and T- {ab} = tp {ab} = (^N {a c N b) d - ^N ab N cd ^j 


4>cd- 


(4.4) 


(4.5) 


In the above, we have introduced the notation T/{ a c i to denote PSTF with respect to n a . As in the 1+3 formalism, 
the existence of the preferred spatial congruence and induced metric N ab allows us to define two derivative operators. 
The first is a convective derivative along the spacelike congruence 




n e D e ip a _ b c 


(4.6) 


and the second is a totally projected derivative on the 2-surfaces 

A := Nj Nj... AV".\V ■■■N i d D i ip f g h - i . (4.7) 

We can now decompose the covariant spatial derivative of n a into its irreducible parts[8, 9] 

Da^b — + ~(pN ab + £ C ab + (4*8) 

where 


o, a — ti D c n a — n a , (4.9) 

(p = 6 a n a , (4.10) 

£ = \e ab S a n b , (4.11) 

C ab = 5{ a n b} . (4.12) 


Here, we interpret <p as the expansion of the 2-surface, (, ab as the distortion of the 2-surface (or shear of n a ), a a the 
acceleration and £ the rotation of n a . Likewise, we can also perform an irreducible split of the timelike derivative of 


n a = Au a + oi a where a a = fi a and A = n a u a . 


(4.13) 


The remaining 1+1+2 variables arise from the irreducible decomposition of the usual kinematical and gravitational 
variables in the 1+3 formalism [8, 9]: 


ii“ = An A + A a (4.14) 

= nn a + n a (4.15) 

o ab = E ^n a n b - ^Nab'j + 2E (a n b ) + E ab (4.16) 

N a b = & {^bia'db ^N ab ^ + 2£^ a n b ^ + £ ab (4.17) 

H a b = 7-L [n a n b — -N ab j + 2H( a n b ) + H a b- (4-18) 
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The key variables in the 1+1+2 formalism for vacuum spacetimes are therefore 

{&,A,(!),X,£,n,H,Z,A a ,n a ,'Z a ,a a ,a a ,£ a ,-H a ,X ab ,£ ab ,-H ab ,U} (4.19) 

In the extension to studying perturbations to the f(R) Schwarzschild black hole, we must also implicitly include the 
1+1+2 energy-momentum variables in order to describe the effective curvature fluid. The 1+1+2 energy-momentum 
tensor can be written as follows: 

T a b f. lU a Ub + p (N a b ^a^b] + 2 Qu^a^b) ~b aQb ) T II i.Tl a Tlb ~^N ab j + Yl^ a Tl b ^ + II ab- (4.20) 


The full decomposition of the curvature fluid energy-momentum tensor into its 1+1+2 form is given by [41]. The 
linearisation procedure around a Schwarzschild proceeds by substituting Eq. (2.11) into Eq. (2.8) and performing a 
1+1+2 decomposition. Zeroth order terms are given by the set of Schwarzschild scalars {£, </>, .4} and all products of 
first order terms or higher are dropped from the resulting expressions, leading to the following definitions [41] 


K = T fZR + Mgii + fZPR 


P R = 


1 - 

TK 

1 

To 


fSR - MUR - \<t>%R - lfoS 2 R - \fo R 


i 


QR = -ri 
Jo L 

Qa=-fSaR, 

J 0 


fo - AR) 


= 


n R = — 


fo 


0 L 


fo 

/(" 


0 L 


\fo&-\<ffoR-\f'T 2 R 


foSaR-^foSaR 


Kb = -jr TMR- 

Jo 


(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 


The identities 


nf Qh } = - S {a Qb } and e c{a n t} c - e c{a 5 c n b} + ^e c{a II fc} c = 0, (4.28) 

have been used to simplify the evolution and propagation equations for the electric and magnetic Weyl 2-tensors as 
given in ??. 


B. The Equations 

The full system of 1+1+2 equations arises from the decomposition of the 1+3 equations coupled with a further 
system of equations defined by the Ricci identities for n a . 


R a bc = 2 V[ a V b jn c — R a bcd,n d = 0. (4.29) 

The Ricci identity applied to n a guarantees that we have sufficient equations in order to determine the 1 + 1 + 2 
variables. In particular, this third rank tensor can be covariantly split with respect to u a and n a to give dynamical 
equations for the variables derived from the irreducible decomposition of n a (a a , a a , </>, £, £ a b) [8, 9]. The complete 
system of equations can be covariantly split into propagation, evolution and constraint equations by projections with 
respect to u a , n a , N ab as well as various symmetry operations. Evolution equations involve dot derivatives, propagation 
equations involve hat derivatives, mixed propagation and evolution involve both and constraint equations are restricted 
to the 2-surfaces. The full system of 1+1+2 equations can be found in [8, 9] and the system of 1+1+2 equations for 
the f(R) Schwarzschild spacetime in [41]. We will not reproduce these equations here. 












V. The Regge-Wheeler Tensor 



A. Review: The Background Spacetime 


The Schwarzschild spacetime is covariantly characterised by the following non-zero 1+1+2 scalars: {A, £,</)}. The 
system of equations that covariantly describes the Schwarzschild spacetime can be written in the following compact 
form: 


$=-l<j?-£, 

(5.1) 


£ = —£(/>, 

(5.2) 


A — —A ( (f> + A '), 

(5.3) 

with the constraint equation 

£ + cj)A = 0. 

(5.4) 


A parametric solution for these equations can be found [9]: 


/ 2 M 

„ M 

2 M1 

\/l-, 

A= 

1 - 

V r 


r 


£ = - 


2 M 


(5.5) 


B. Review: Gauge Invariant Perturbations 

In more conventional approaches to linear perturbation theory in GR, the perturbation to some tensor quantity 
is the difference between its value at some event in the physical spacetime and its value at the corresponding event 
in the background spacetime. Following [49], if we consider a background 4-dimensional spacetime Mo then we can 
introduce a 5-dimensional smooth manifold that contains a 1-parameter family of smooth non-intersecting 4-manifolds 
M e - These manifolds are interpreted as perturbations to the background spacetime. Implicitly, this means that to 
define a perturbation we must choose a gauge, or a mapping from Mo to Me- If Q is some geometric object in the 
background and Q e the corresponding object in the perturbed spacetime then we can define a Taylor series expansion 
along the integral curves of a vector field V that is everywhere orthogonal to M e via the 1-parameter group of 
transformations h e from Mq to M e - 


Qe — he* [Qo + e {C-vQe) e =0 + ^(e 2 )] ■ (5-6) 

The immediate implication here is that the first order term, the linearised perturbation of Q, will be gauge-invariant 
iff 


C-vQ o = 0, 


(5.7) 


for all 4-vectors V on Mq. This will be the case if: 

• Qo vanishes in the background. 

• Qo is a constant scalar field. 

• Qo is some linear combination of Kronecker deltas with constant coefficients. 

This is known as the Stewart-Walker lemma [48] and provides a very powerful foundation from which gauge-invariant 
variables can be defined in the covariant approach to perturbation theory. This is due to the fact that in the covariant 
approach we aim to avoid any reference to a background spacetime and only invoke the background when determining 
which variables are zeroth order. Gauge-invariance of our set of equations is then guaranteed by the Stewart-Walker 
lemma if we find a complete set of variables that vanishes in the background. For the Schwarzschild spacetime, the 
only non-zero variables in the background were {£,</>, A} with all other scalars, vectors and tensors vanishing in the 
background due to symmetry. To construct gauge-invariant variables from these non-vanishing scalars, we just need 
to invoke spherical symmetry, which necessitates that the angular derivatives of the background scalars must vanish 
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and will therefore constitute a set of gauge invariant variables. The set of variables which we will therefore work with 
in the linearised perturbation equations are as follows: 

5 a {£,4>,A} = {X a ,Y a ,Z a }. (5.8) 

In the formalism we have outlined above, we also have another implicit type of ’gauge 1 freedom, namely the freedom 
to choose a frame basis in the tangent space at each point. This is the same as in any other tetrad based approach, 
such as the Newman-Penrose formalism. As in [9], we will use the term gauge-invariance to refer to invariance under 
a mapping between the true and background spacetimes and frame invariance to refer to invariance under our choice 
of frame vectors. 


C. Review: The Regge Wheeler Tensor 

In the metric approach to GR, a master equation for axial gravitational perturbations was originally derived by 
Regge and Wheeler [45] with the corresponding polar equation being derived by Zerilli [56] over a decade later. The 
structure of these wave equations can be shown to reduce to a Schrodinger like equation defined in terms of an 
effective potential V e s. These original studies have formed the basis for many different approaches to studying linear 
perturbations to the Schwarzschild black hole [9, 10, 31, 32, 34, 36, 44, 54]. In the 1+1+2 approach to gravitational 
perturbations of a Schwarzschild black hole, the master variables for the polar and axial sectors can be unified into a 
single covariant, gauge invariant and frame invariant tensor Wy ab y [9]. 

Building on this work, a master variable Wy ab y for gravitational perturbations for the f(R) Schwarzschild black 
hole has recently been derived in [41] and is given in terms of the shear 2-tensor (/ ab , a TT tensor constructed from 
the radial part of the electric Weyl tensor Sy a X b \ and a TT tensor constructed from the curvature scalar 5y a S b yR. 
The dimensionless, gauge-invariant, frame-invariant and transverse-traceless tensor W ab is defined by [41] 1 

W ab = ^r 2 C, ab - + ^r 2 S {a 6 b} R, (5.9) 

where the even part of W ab is coupled to a curvature term: 

\ f jrJ { a5 b} R =\^ b . (5.10) 

Consequentially, we have to include the trace equation for the curvature scalar in order for the wave equation to close. 
The curvature term vanishes in the axial sector and the variable W ab reduces to exactly the same form as in GR. The 
covariant wave equation this variable obeys is given by [41] 

W {ab} - & {ab} - AW {ab} + (, i> 2 - £) W {ab} - 6 2 W {ab} = 0, (5.11) 

which is identical to the covariant wave equation in GR [9, 43]. Note that both the even and odd parity components 
of W ab obey the same wave equation. This tensor equation can be decomposed into scalar harmonics (see Appendix 

C) 


W-W-AW 


'e(£ + 1 ) 

r 2 


+ 3£ 


W = 0, 


(5.12) 


where W = {Wt,Wt}- Associating the hat derivative with an affine parameter p. i.e. ~ = d/dp , we can convert to 
the parameter r, p —> r, and then switch to the tortoise coordinates of Regge and Wheeler 


r,=r + 2Mln(^-l). 


(5.13) 


Letting 


= V’rw = w , 


(5.14) 


1 Note that [41] denotes this tensor as M a & whereas we use the notation of [9]. 
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we find that the harmonic equation reduces to a Schrodinger type equation 


with an effective potential Vr 


Vt 



e(e + i) 

r 2 


6M 

r 3 ’ 


(5.15) 


(5.16) 


which is just the Regge-Wheeler (RW) potential for gravitational perturbations to the Schwarzschild spacetime [9, 
41, 45]. 


D. Review: Scalar Perturbations 


In FOG we note the presence of scalar modes that are not possible in Einstein GR. The trace equation yields 
a wave equation in terms of the Ricci scalar R. The equation obeys the same functional form as the generalised 
Regge-Wheeler equation for massive scalar perturbations in Einstein GR but with an effective mass given by [41] 


U 2 


A 

3 fs 


The covariant wave equation that the variable obeys is given by [41] 


it - n - An - (e - u 2 + s 2 ) n = o, 

where 7 Z = r R. Introducing scalar spherical harmonics, this equation reduces to [41] 


A ~ TZ-S — Alls 


£-U 2 


t(t + iy 

r 2 


TZs = 0 , 


(5.17) 


(5.18) 


(5.19) 


where U = Ci/^C^), with C\ and C 2 constants. Converting from the parameter p to r and then introducing tortoise 
coordinates, we find [41] 


where 





n = 0, 




^+ 1 ) , 2M f 2 

r 2 ' r 2 ' U 


(5.20) 


(5.21) 


VI. The Weyl Terms 

We now move on to the core subject of the paper, namely the description of gravitational perturbations via the 
Weyl tensor in the 1+1+2 formalism. The definition of a master variable for gravitational perturbations is non-unique 
(see e.g. [36]) and different definitions will provide different advantages or insight. The previous definition of the RW 
tensor W a b given in Eq. (5.9) allows us to interpret the RW equation as being related to the fluctuations in the radial 
part of the electric Weyl tensor as well as the distortions of the sheet, given by ^ a f,. Alternatively, we may directly 
invoke the decomposition of the Weyl tensor to introduce a set of new master variables. This allows us to relate 
the 1+1+2 RW tensor to h + and h x , the NP Weyl scalars and the 2+2 gauge invariant master variables in a clear, 
transparent manner. 


1. Unified Polar and Axial Gravitational Perturbations 

It is a well established fact that gravitational waves propagate in a perturbed Schwarzschild spacetime and it would 
be reasonable to expect that the transverse-traceless tensors governing the gravitational perturbations will obey wave 
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equations that reduce to the plane wave case in the limit M —>• 0. The most natural variables to study would be 
those relating to the transverse-traceless (TT) degrees of freedom in the Weyl curvature tensor, as this provides a 
description of the free gravitational field. In the 1+1+2 formalism, the TT components of the electric and magnetic 
Weyl tensor will just be the 2-tensors £ a b and H a b- Constructing the wave equations for these variables, however, we 
find that they do not form a closed system, they contain forcing terms from other 1+1+2 tensors. Such couplings 
are not present in the plane wave limit making the interpretation and solution of these wave equations in isolation 
non-trivial. 

In order to find a wave equation for the Weyl curvature 2-tensors, we need to perform a 1+1+2 decomposition 
of the 1+3 evolution equations for the Weyl curvature (see Eqns. 1.3.42 and 1.3.43 in [51]). Linearising about the 
Schwarzschild background, this yields [9, 41]: 

£{ab} - e C {a^ b } C = -<+{a<5 C ?4} + e c{a?4} C ~ ^^ab, (6.1) 

7~t{ab} + e c{a£ b y C = e c{a^°^b} + ^^{aCb} ~ ^2^ ^ tc ^ b } C ' (6-2) 


The covariant wave equations for the electric and magnetic Weyl 2-tensors can now be derived by applying the wave 
operator 4' — 'L to £ a b and Hab using the above expressions to yield 


R{ab} - ?t£{ ab } — (0 + 5-4) H{ab} ~ 5 2 H{ab } + 
£{ab} £{ab} (0 T 5-4) £{ab} ^ £{ab} T 






n [ab} = [-3£ e {a £ b} L + 25 {a H b} \ (0 - 2-4) (6.3) 

£{ab} = [3£C{ab} + 2 5 {a £ b} ] (0 - 2-4). (6.4) 


Clearly, these wave equations are coupled to the 1 + 1+2 tensors £{ a f>}> C{ab}, £a and H a • Neither of these variables 
alone will obey a closed system of equations. By inspection of the propagation and evolution equation we can start 
to see similarities in the structures of these two wave equations. Importantly, the coupling terms can be on the RHS 
of Eqns. (6.1) and (6.2) share similar pre-factors that are reflected in the RHS of Eqns. (6.3) and (6.4) suggesting 
that we should be able to decouple these 2 -tensors from the other 1 + 1+2 variables. 

Following this insight, we can construct a new variable from the linear combination of the electric and magnetic 
Weyl 2-tensors J± h that completely decouples from the other 1+1+2 variables. Explicitly, this new perturbation 
variable is given by 


3{ab} = £ {ab} ± e C {a?4 } C , (6-5) 

= \£t + Ht\ Qab + \£t ± Ht\ Qab, (6.6) 

where in the second line we have decomposed the perturbation variable into tensorial harmonics as per Appendix 
C. More importantly, this variable can be shown to obey the following closed, covariant and gauge invariant wave 
equation 


^6} - AU - (•* + W Jtab} T (4-4 - 24>) j± 6} 

which, when decomposed into covariant harmonics, reduces to 
J ± - i ± - (-4 + 30) J ± + (4-4 - 20) j ± + 


-{5 2 + 2K)J± b} + ^A 2 -±£)jf ab} 


30 2 + 2 £+ 4 - 4 2 

r z 2 


= 0 , 


= 0 , 


(6.7) 

( 6 . 8 ) 


where = {-7+, }. As with W{ ab }, this variable unifies the even and odd parity perturbations into a single 

transverse-traceless tensor J ab . By construction, this variable is defined by the transverse-traceless degrees of freedom 
in the electric and magnetic Weyl tensors. It is therefore natural that this variable should describe the propagation of 
gravitational waves in the Schwarzschild spacetime. In the limit M —> 0, we recover the plane wave case, as expected. 
Additionally, we note that gauge invariance of our perturbation variable is guaranteed by the Stewart-Walker lemma 
[48] as all vectors and tensors vanish in the background spacetime. 

In the context of f(R) gravity, it may seem a little curious that the wave equation for this perturbation variable 
is exactly the same as found in General Relativity. However, this result should not come as too much of a surprise. 
The perturbation variable, being TT in nature, explicitly describes massless modes that propagate along null curves. 
The scalar modes, having an effective mass of U 2 , will physically propagate along timelike curves on the black 
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hole background. In addition, we are explicitly considering a constrained class of solutions in f(R) for which the 
scalar modes are not excited in the background spacetime. As such, the scalar modes and tensor modes decouple at 
linear order. This decoupling between the massive and massless modes has previously been noted [35, 41] but the 
derivation via metric based approaches can often obfuscate the underlying physics. By adopting the geometrically 
and physically meaningful 1+1+2 semi-tetrad approach, we can explicitly see which kinematical and gravitational 
terms are of importance and we can isolate physically meaningful variables, such as the TT components of the electric 
and magnetic Weyl tensors. 

As a side note, by geometrically and physically meaningful we are referring to the notion that everything can be 
related to the dynamics and kinematics of congruences as well as the fact that all objects are gauge-invariant and 
covariantly defined without reference to a given coordinate system. The derivative operators appearing in the wave 
equation are just convective derivatives along the timelike and spacelike congruences or angular derivatives on the 
2-surface. The perturbation variable itself tells us about the evolution of the Weyl 2-tensors that are defined on the 
2-surfaces and therefore provides a natural description of gravitational radiation. In addition, these variables can be 
transparently related to the broader class of LRS-II spacetimes, where it can be shown that the perturbation variable 
defined above describes gravitational perturbations to vacuum LRS-II spacetimes [? ]. 

Finally, it is also worth noting that were we to consider more astrophysical realistic systems or systems with 
matter content, the role of the scalar modes can become pronounced. An example of such behaviour can be seen 
in the dynamical secularisation of isolated neutron stars in the late-inspiral and merger phases of a compact binary 
coalescence [1]. A study of neutron star systems in the 1+1+2 formalism would provide a natural extension to this 
work. 


2. Axial Gravitational Perturbations 

As mentioned previously, the curvature term only occurs in the even parity sector. This means that axial gravita¬ 
tional perturbations to the f(R) Schwarzschild black hole will be governed by the same covariant wave equation as in 
General Relativity. A natural candidate for such a perturbation variable is the magnetic Weyl tensor. Consequentially, 
we introduce a perturbation variable V{ a b} constructed from the radial part of the magnetic Weyl tensor 

V {ab} =r 2 6 {a S b} n. (6.9) 

This variable can be shown to obey the following covariant wave equation 


^{ab} — V{ab} — (A + 3(/>) V{ab} ~ + 2A ] V{ a f,} — 0. 

Decomposing the above equation into tensor harmonics 2 we find that 


V - V - (A + 3^) V + 


e(£ + i) 3 2 ■ 

A 2 2*^ +6£ 


v = o. 


( 6 . 10 ) 


( 6 . 11 ) 


Though this harmonic equation looks different to the Regge-Wheeler equation, this variable is related to Regge- 
Wheeler variable. This can be seen by rescaling the tensor defined above as follows: X = r 3 V, Eq. (6.11) then reduces 
to 


0 = 


X - X - AX + 


l{l + 1 ) 


38 \X 


= X - X - AX 


£{e + i) 

r.2 


3 e\x. 


( 6 . 12 ) 

(6.13) 


This is nothing more than the Regge-Wheeler equation as per Eq. (5.12). This allows us to reinterpret the RW master 
variable as simply the radial part of the magnetic Weyl scalar as R = n a n b H a i,. This is more physically intuitive than 
the original derivation via metric perturbation theory and is manifestly covariant and gauge invariant. 


2 As per Appendix C 










13 


VII. Relation to Other Formalisms 
A. The Newman-Penrose Formalism 


Another way to interpret our new perturbation variables is to re-express the Newman-Penrose scalars in terms of 
the 1+1+2 variables. The Newman-Penrose formalism is a full tetrad approach in which the underlying frame vectors 
form a null tetrad consisting of two real null vectors ( l a ,k a ) and a complex-conjugate pair ( m a ,fh a ) [37]. The two 
real null vectors correspond to ingoing and outgoing null congruences, whereas the complex-conjugate pair will simply 
correspond to a decomposition of the 2-surface orthogonal to u a and n a . The frame vectors are defined as follows: 


l 

k 


a 


a 


m 


a 


l 

71 

i 

71 

i 

71 


( U a + n a ) 

(Ua - n a ) 

(v a - iw a ) 


with the metric being decomposed as follows 


U a = 0 

k a k a = 0 l a k a = -1 

m a rh a = 1 m a m a = rh a fh a = 0 , 


9ab — lak b k a l b + 

= ~u a u b + n a n b + N ab , 


(7.1) 

(7.2) 

(7.3) 

(7.4) 

(7.5) 


where we have identified the induced 2-metric with the two complex frame vectors: N ab = 2m( a fh b y The Newman- 
Penrose scalars are defined by the appropriate contractions of the Weyl tensor with respect to this null tetrad and 
are given by 3 


= C abcd l a m b l c m d (7.6) 

* i = C abcd l a m b l c k d (7.7) 

4 > 2 = C abcd l a m b fh c k d (7.8) 

*3 = C abcd l a k b fh c k d (7.9) 

*4 = C abcd m a k b m c k d , (7.10) 


with the 5 complex Weyl scalars encoding the 10 independent components of the Weyl tensor. The conventional 
physical intuition applied to these scalars is as follows: 4+ describes transverse radiation along k a and thereby 
describes ingoing gravitational radiation. 4*i is an ingoing gauge wave described by the longitudinal radiation along 
k a . 4*2 is a Coulomb term related to gravitation attraction and frame dragging effects. +3 describes longitudinal 
radiation along l a and is therefore an outgoing gauge wave. Finally, 4>4 describes transverse radiation along l a and is 
thereby characterises outgoing gravitational radiation. 

One of the motivations for considering the Weyl scalars is that the response of a GW detector will be encoded in 
these curvature scalars in the Jordan frame [16]. In particular, for outgoing waves, the Weyl peeling theorem tells us 
how the Weyl scalars fall off along outgoing radial null geodesics in some neighbourhood of future null infinity J + : 

y n ~r n ~ 5 . (7.11) 


This result is relatively robust and valid for a wide range of tetrad frames given a suitable choice for r. In an 
asymptotically flat spacetime, given that we are sufficiently far away from the source, the outgoing waves can be 
expressed within the realms of linearised gravity in the TT gauge: 

4-o = 4b = 4> 2 = ^3 = 0 4-4 = -h + +ih x , (7.12) 


with h+ and h x the standard gravitational wave polarisations, i.e the 2 graviton degrees of freedom. This means that 
in General Relativity, the radiative degrees of freedom that decay as r _1 and are observable by GW detectors are 
simply encoded in 4 , 4 . The extension to f(R) gives us the possibility of exciting an additional transverse, radiative 
scalar mode. However, under many general considerations, the radiative component of this scalar mode tends to 


3 Note that other definitions or conventions exist. 
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vanish due to constraints on the underlying scalar-tensor theory [ 1 ]. As such, the scalar modes will only couple 
weakly to a GW detector making direct detections difficult [1, 13]. 

Performing a systematic decomposition of all the Weyl tensor into it’s constituent 1+1+2 components, the Weyl 
scalars can be re-written as follows [47]: 


4+ 
'hi 
'h 2 
4fo 


T 4 


£ab 4“ G-f a 77 




m a m b 


--^= [£ a - <iabH h ] n 

\\£-m 

[£a + tabH b ] ffi a 


£ab £r{a/H' 


b} 


m a m b . 


(7.13) 

(7.14) 

(7.15) 

(7.16) 

(7.17) 


This decomposition follows from analogous results derived in a 3+1 or 1+3 approach [2, 3, 47]. As is usual in 
a tetrad formalism, these scalars will, in general, be frame-dependent with the behaviour of these scalars under 
transformations explicitly known [49]. Reassuringly, the only Newman-Penrose scalar that is non-vanishing in the 
background spacetime is 4 '2 for which we recover: 



M 

+3 ’ 


(7.18) 


as per the parametric solution for the 1+1+2 variables in Eq. (5.5). The above results are also consistent with the 
literature in the sense that all LRS spacetimes will be of Petrov type D or O [17] 4 . This can be seen by the simple 
requirement that under LRS, all 2-vectors and 2-tensors must vanish in the background. This leaves 4^ 7 ^ 0 with the 
Petrov type determined by £ and 77. If these scalars vanish the Petrov type is O, otherwise the Petrov type is D. 
Under perturbations, the Petrov type will be asymptotically type N due to the Weyl peeling theorem. We can extend 
this discussion to include the four non-zero vacuum invariants constructed from the Riemann curvature tensor [2] 

1 = 34^ -44'i4' 3 + 4+4^, (7.19) 

J = 4' 0 4' 2 4' 4 + 24'i 4' 2 4' 3 - 4'o4’g - 4' 2 4' 4 - 4>2- (7.20) 

In a Petrov type D spacetime, as Petrov type O is trivial, these invariants reduce to 

I = 34 , 2 and J = -4^. (7.21) 

From our decomposition of the NP scalars, these can be expressed via the 1+1+2 scalars as 


I = l [£■ 2 - n 2 - 2i£H] , (7.22) 

J = -1 [£■ 2 - 3i£ 2 n - 3£77 2 + iH 3 ] , (7.23) 

8 

J 2 

S = 27—r = 1 For Type D or Type II Spacetimes. (7.24) 


From these decompositions we can immediately identify the correspondence between the perturbation variable V a b 
and the imaginary part of 4 , 2 : 

S[r 3 4' 2 ] ~ ipRw ~ r 3 77. (7.25) 


The relationship between the imaginary part of 4t 2 and the Regge-Wheeler variable was first discussed by [44] 5 . 

Similarly, we can immediately identify a correspondence between {4' 0 , 4' 4 } and the perturbation variable J^ h 
Notably, we see that 


^0 = JX m a m b , (7.26) 

^4 = J~ b m a fh b . (7.27) 


4 See Appendix A for a brief review of the Petrov classification. 

5 Note that [44] uses a different notation for the NP scalars based on their spin-weight: {4 f _2,vE f — i,4/o,4/i,4/2}- 
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This allows us to associate 3A with outgoing gravitational radiation 6 and reinforces the notion that 3A, contains, in 
a compact form, the gravitational wave polarisations {h + , h x }- 

By construction, the scalar degrees of freedom present in the class of f(R) theories considered here do not enter the 
Weyl scalars. This is reflected in our ability to both construct a closed, covariant wave equation that decouples from 
the scalar modes as well as the inability of the scalar modes to affect the potential of the tensor degrees of freedom 
at linear order. 

A useful approach for studying the curvature-fluid perturbations is via the Ricci scalars. In the most general case, 
these scalar quantities will encode the ten independent components of the Ricci tensor in the form of three real scalars 
{$ 00 ! $ 11 , $ 22 } and six complex scalars {$01 = $io, $02 = $ 20 , $12 = $ 21 }- Unlike the Weyl scalars, the Ricci scalars 
are implicitly related to the energy-momentum distribution of the spacetime via EFE. In the context of perturbations 
to the f(R) Schwarzschild black hole, these scalars will measure excitations of the curvature fluid. The Ricci scalars 
are defined as follows [47]: 


\R ab l a l b i 

(7.28) 

\R ab ( l a k b + m a fh b ), 

(7.29) 

\Rabk a k b , 

(7.30) 

$10 = - \R ab l a m b , 

(7.31) 

$20 = ^ R ab m a m b 

(7.32) 

$21 = \R ab fh a k b . 

(7.33) 


As per the Weyl scalars, we can systematically decompose the Ricci scalars into their constituent 1+1+2 variables. 
We find the following: 


$00 

$11 

$22 

$01 

$02 

$12 


\{p-Q + \ n ), 

(7.34) 

^ (n - n + n ab m a m b ) , 

(7.35) 

\ (p+ Q + l n ). 

(7.36) 

2V2 ( ' IIa ^ m ’ 

(7.37) 

^11 ab m a m b , 

(7.38) 

2^ *5“) m ■ 

(7.39) 


The 1+1+2 expressions presented above are valid in both General Relativity as well as the f(R) extension considered in 
this paper. In the Schwarzschild background, these variables naturally vanish as we only consider vacuum gravitational 
perturbations, i.e. T^> vanishes to all orders. In the f(R) extension however, these variables will vanish in the 
background but will be non-zero once we consider first-order perturbations involving the scalar degrees of freedom. 
For example, we can explicitly insert the linearised energy-momentum tensor for the curvature fluid into $00 and $22 
to obtain: 


fo 

R - \r - R 


- h 2 R + R~- 

AR 

V'o 

3 

K 6 / 

6 


fo 

R+-R-R 

^a + -A 

- -S 2 R--A 


Vo 

3 

K 6 / 

6 3 



(7.40) 

(7.41) 


6 See Appendix B. 
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B. The 2+2 Formalism 


The other approach to which we wish to make a connection is that of the 2+2 formalism. This formalism was 
originally introduced in [20, 21] and a completion of this formalism was given by [22] who providing a systematic 
derivation of the gauge-invariants, introduced a fluid-frame decomposition and wrote down the resulting system of 
master equations governing gravitational perturbations of spherically symmetric spacetimes. The approach is built 
around the decomposition of the background 4-dimensional spacetime Ai 4 into a warped product Ai 4 = Ai 2 0 S 2 , 
where Ai 2 is a 2-dinrensional Lorentzian manifold and S 2 is the 2-sphere. In essence, we are covariantly factorising out 
the spherical symmetry and reducing the problem to a two-dimensional problem written in terms of the two essential 
coordinates: time and radius. As such, the metric can be written as' 

ds 2 = gAB{x C )dx a dx B + r 2 (x c ) r y a bdx a dx b , (7.42) 


where gAB is the metric on Ai 2 and j a b is the metric on the unit sphere S 2 . The scalar r = r(X a ) is defined on 
Ai 2 and can be identified as the invariantly defined radial coordinate of spherically-symmetric spacetimes. Using this 
decomposition, the EFE in vacua can be re-expressed as 


Gab = —2 (va\b + vavb) + 2 + 3 vcv c - gAB = tAB, 

\g\ = vJ c + v c v c -n = Q, 


(7.43) 

(7.44) 


where 1Z = \R A is the Gaussian curvature of Ai 2 . The energy-momentum conservation equations take the following 
form: 

t AB lB + 2 t AB v B - 2 v a Q = 0. (7.45) 

Following [22], all metric perturbations can be written as a scalar, vector or tensor field on Ai 2 times a spherical 
harmonic scalar, vector or tensor field on S 2 . As such, the axial metric perturbations can be written as: 


h Axial = (0 ^Axial \ 

^ hY ab ) 


(7.46) 


and the polar metric perturbations as: 


t Polar 


h AB Y /jPoiar Y a 

^Polar^ r 2 (KY lab + GY ab ) 


(7.47) 


Note that there is a natural correspondence between the covariant harmonics {Q, Q a , Q ab } and their spherical harmonic 
counterparts {Y, Y a , Y ab }, this is discussed in Appendix C. In practice it is often convenient to adopt the Regge-Wheeler 
gauge {h, h^ lar , G} = 0 as, in this gauge, the remaining gauge-invariant variables simply correspond one-to-one with 
the bare perturbations. Knowing how the variables transform under general gauge-transformations, we could always 
map to a different gauge using known algebraic relations. 


1. Polar Perturbations and Correspondence 

Using a radial-frame in Ai 2 , the metric perturbation kAB can be split into three gauge-invariant scalars {\, <£>,?} 

kAB = (x + P) ( n A n B + u A u B ) + ? («a«b + n A u B ) ■ (7.48) 

Substituting this decomposition into the 2+2 EFE, we can recover a system of scalar equations by the appropriate 
projection operations with respect to the basis { u , n}. The system of master equations for the Schwarzschild spacetime 


7 Note carefully: In this section, {A, B,... } denote coordinates on M 2 and (a, b,... } denote coordinates on S 2 . 
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in GR can be written as follows [22, 36] 8 


^X + X = -2 


2v z - 


6 M 


(X + <P)~ - 2W) x + 


/ , + 2 )(£ — 1 ) 


, ,4m. , (* + 2)(*-l) 

= -WX - ^ - 7T (X + ¥>)- ^2 - Z X, 

-? = 2^ (x + +) + x'- 

In addition, we also have three constraint equations on the Cauchy data for {x, </>,?}: 

=-wx+ e ^ 2 1 \ , 

<p" = (X + ¥>) - (g + 2 2 ) r ( 2 ~ 1) X + Wx' - 2W<f/, 


s’ = ~2vs - X - 2 ip. 


(7.49) 


(7.50) 

(7.51) 

(7.52) 


Clearly, from the above, the highest derivatives of x form a wave equation with characteristics set by the metric g Mi , 
and Cauchy data {x,x} that can be set independently of matter perturbations [22]. It was therefore noted in [22] 
that x can be reasonably said to characterise polar gravitational waves. The interpretation of <p is a less trivial but 
was originally noted to correspond to longitudinal gravitational waves made physical by the presence of matter. As 
we work in a vacuum spacetime, how we should interpret ip is not clear. Finally, <; was noted to be a term that 
is advected with the fluid. These physical descriptions of the gauge-invariants can be made manifest by explicitly 
writing down a correspondence between the 1+1+2 variables and the 2+2 gauge invariant perturbations. In essence, 
this approach to understanding the physical and geometrical meaning of gauge invariant perturbations is reminiscent 
of [ 6 ] but adapted to spherically symmetric spacetimes. 

Performing a systematic study of the 1+1+2 variables in terms of the 2+2 gauge invariants, we find the that the 
gauge-invariant perturbations are captured by the electric and magnetic Weyl 2-tensors: 

£a b = ~l(x+<p)Y a 6 , (7.53) 

n ab = -^Y ab . (7.54) 

In many ways this is highly reassuring. After all, it was noted that x and <p somehow relate to polar gravitational 
wave degrees of freedom. Similarly, the interpretation of g was less transparent but, via the 1+1+2 approach, we can 
clearly see that this gauge-invariant perturbation is related to the magnetic Weyl 2-tensor and therefore will be related 
to genuine relativistic effects, such as frame-dragging. This also corresponds to the the result presented earlier in 
Eq. (6.7) whereby the TT degrees of freedom of the Weyl tensor were shown to obey a closed covariant wave equation. 
As such, the master variable an d its concomitant wave equation are simply a covariant, gauge-invariant and 

frame-invariant repackaging of the system of equations given by Eq. (7.49) for the gauge-invariant perturbations 

{x,++}- 

From the result in Eq. (6.7), which is valid for both General Relativity as well as /(R), it would be natural to 
assume that were a 2+2 decomposition of the f(R) field equations to be performed, an analogous closed system of 
equations for these gauge invariant perturbations could be found. Implicitly, we would also have to fold in first-order 
energy momentum perturbations corresponding to excitations of the scalar degrees of freedom. This would follow a 
procedure analogous to that used in the 1 + 1+2 decomposition and linearisation of the energy-momentum tensor for 
the curvature fluid. 

As a final point, we can complete the correspondence between the 1+1+2, 2+2 and NP polar sector perturbations 
by studying the perturbations to fI /4 in terms of the 2+2 gauge-invariant variables. In the background spacetime, ^4 


Note that in this section a dot derivative is defined by u a Va and a prime derivative by n A S7 a, where Va is the covariant derivative 
on M 2 , i.e. V C 9AS = 0. 
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vanishes and therefore its perturbation is naturally gauge-invariant. This perturbed Weyl scalar can be written as 


5^4 = 5C abcd fh a k b fh c k d 

(7.55) 

= 5RaBcDffi a fh C k B k D 

(7.56) 

= -^k B k D k BD [rh a m c Y ac ] 

(7.57) 

= ~ \ bi + <P + 4] [m a fh c Y ac ] 

(7.58) 


(7.59) 


The expression in the last line is obtained by decomposing the radial-frame metric perturbation k AB into the 2+2 
gauge-invariants as per Eq. (7.48). Immediately we can see that this expression is nothing more than the harmonic 
decomposition of our perturbation variable J~ c projected with ffi a fh c . This is in agreement with the correspondence 
between the 1+1+2 master variables and the NP Weyl scalars as well as the correspondence between the 1+1+2 
variables and the 2+2 gauge invariant perturbations. 


2. Axial Perturbations and Correspondence 


In exactly the same manner as we did for the polar perturbations, we can analyse the axial sector and how the 
results correspond to those derived in both Eq. (6.7) and Eq. (6.10). Luckily, the analysis is greatly simplified as the 
system of equations for the axial sector are less entangled than their polar sector counterparts. Perturbations to the 
axial sector, in vacuum, are completely characterised by a single scalar variable II defined by [20-22] 


n = 


AB 


'k A ) 


- 1 - 


I B 
2 M 


- 1/2 


k' t —kt + 
r 


1_2MX 


1/2 


(7.60) 

(7.61) 


where we have explicitly evaluated II in terms of its metric components for Schwarzschild coordinates. The corre¬ 
sponding scalar wave equation is 


1 

r 2 



(£ + 2 ) (£-l) = 0 . 


(7.62) 


On a Cauchy surface, the two first-order degrees of freedom are given by {II, n}. Explicitly evaluating the covariant 
derivatives, we would find a scalar wave equation with characteristics set by the metric: 


-ii + n" = s n 


(7.63) 


Consequentially, it is reasonable to state that II describes axial gravitational waves. Expressing the 1+1+2 variables 
in terms of 2+2 gauge-invariants we can immediately identify the axial master variable with the radial part of the 
magnetic Weyl 2-tensor 


U = 


1(1 +l) n 
2 


(7.64) 


Following our previous discussion, the master variable V{ a b} is a covariant, gauge-invariant repackaging of the 2+2 
gauge invariant master variable II whose physical significance is made manifest by the identification of V a b with the 
radial part of the magnetic Weyl 2-tensor. 

Finally, we can complete the correspondence by noting the link between the perturbed Weyl scalars and the 2+2 
gauge invariants [38], the most significant of which is the perturbation to ^ 2 : 


<54-2 = -^ + i)ny 

~ +RW- 


(7.65) 

(7.66) 

(7.67) 


This summarises the results detailed in this paper, namely that the radial part of the magnetic Weyl Tensor corresponds 
to the Regge-Wheeler variable and this can be expressed in many different ways in many different approaches. To 
make life more difficult, it is clear from Eq. (6.6)) and Eq. (6.9)) that there is a non-uniqueness in the definition of 
a master variable. As the scalar mode R vanishes in the axial sector, these equations and correspondences will hold 
even in the extension to f(R) gravity. 
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VIII. Summary 

We have presented a further analysis of linear, non-spherical perturbations to the Schwarzschild black hole in 
f(R) gravity. Notably, we have advocated using the Weyl tensor as a powerful object for studying the evolution of 
gravitational waves. The two main perturbation variables that we introduce, J{ a b} and V{ a bj, obey closed, covariant, 
gauge-invariant and frame-invariant wave equations that happen to be exactly the same in both General Relativity 
as well as to the f(R) extensions considered here. 

The first perturbation variable Js a b\ is constructed from the transverse-traceless degrees of freedom of the Weyl 
tensor and is a particularly convenient way to study the evolution of gravitational perturbations in the Schwarzschild 
spacetime. Imposing a set of basic constraints and imposing that /(0) = 0 and R = 0 in the background, the resulting 
equation decouples from the scalar modes at the linear level. In part this is due to the setup used, namely that we 
consider a vacuum Schwarzschild spacetime for which there are initially no scalar mode excitations. The scalar modes 
in f(R) gravity correspond to massive modes, with the mass of the particles set by the parameters of the theory /'(0) 
and /"(0), that propagate along timelike curves. The pure tensor modes described by J a b, however, will be massless 
and propagate along null curves. Were we to consider a more complicated spacetime, by either reducing symmetries 
or by introducing a matter to the system, it is likely that the scalar modes and tensor modes will be more entangled. 
As such, our results are in agreement with the literature on the subject [15, 35, 41]. The variable J{ a b} was shown to 
correspond to the radiative degrees of freedom in the Newman-Penrose formalism [37], i.e. \&4 and T 0 , as well as the 
gauge-invariant master variables \ and + in the 2+2 formalism [22]. This should not be of too much surprise, given 
that the NP Weyl scalars are projections of the Weyl tensor, which describes the free gravitational field and one of 
the primary reasons as to why we advocate its use in the 1+1+2 formalism. 

The second perturbation variable V{ a b} is defined to be the radial part of the magnetic Weyl tensor and is a purely 
axial variable. Given that the scalar modes are even parity in nature, it is unsurprising that the form for this equation 
is exactly the same as in General Relativity. In fact, by an appropriate rescaling, it can be shown that the radial part 
of the magnetic Weyl tensor is exactly the Regge-Wheeler variable. Alternatively, it can be shown that this variable 
is exactly 4/2 in the Newman-Penrose formalism [37] or II in the 2+2 decomposition [22]. 

Finally, we highlighted how the curvature-fluid terms can be expressed via the NP Ricci scalars, e.g. $00 and <f>n, 
which would necessarily be zero in General Relativity. 
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A. Petrov Classification 


The Petrov classification allows us to invariantly divide the gravitational fields into six distinct types. The two main 
approaches to this problem, which are completely equivalent, are as an eigenvalue problem applied to the Weyl tensor 
or by looking at the principal null directions. The eigenvalue approach is tantamount to finding a set of eigenvalues 
A and eigenbivectors X ab such that [47] 


l -C abcd X cd = \X ab . (Al) 

The eigenbivectors correspond to the principal null directions and, due to the inherent symmetries of the Weyl tensor, 
will belong to a four-dimensional space of antisymmetric bivectors. Using a null tetrad as per Newman-Penrose 
formalism, we can introduce the following set of bivectors 

Uab — ^[a^b]) 

V a b k[ a Tfl b ] i 
W ab = m [a m b] - k[J b] . 


(A2) 

(A3) 

(A4) 
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Type I 



Type III-* Type N-» Type O 


FIG. 1. Conventional diagram showing the degenerations of the Petrov type [47, 49]. 


This allows us to express the Weyl tensor in terms of the NP Weyl scalars T n 

C abcd = T 0 (U ab U cd ) + Tr ( U ab W cd + W ab U cd ) + T 2 (V ab U cd + U ab V cd + W ab W cd ) (A5) 

+ 4-3 (V ab W cd + W ab V cd ) + T 4 (V ab V cd ). 


The Petrov classification can now be distinctly written in terms of which NP Weyl scalars vanish. The classification 
is as follows: 
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In terms of the Petrov type, the Weyl peeling theorem can be re-expressed as follows 


(A6) 


It should be noted, however, that the peeling behaviour only holds for spacetimes that are weakly asymptotically 
simple and can be extended smoothly to ^ , though this may be too strict in certain cases [28]. 


B. Radiated Energy 


The correspondence between the electric and magnetic Weyl 2-tensors and the NP Weyl scalar T 4 allows us to 
look at the normal equations for the radiated energy flux and radiated momenta. For a 2-surface at infinity S 2 , the 
radiated energy flux is related to T 4 as follows 


dE 

—— = lim 

at r— >00 


= lim 

r—t 00 



(Bl) 

(B2) 


Likewise, the radiated linear momenta can be expressed via a radial unit vector r* = (sin 9 cos (f>, sin 9 sine/), cos 6) in 
flat space as follows 


= lim I —— [ dtlf, 
dt r->oo I 167 T J s 2 


= lim < - / dfl fi 

r-too I 167 T J S 2 



(B3) 


rt 2 1 

/ J{a b }™ a rn b dt' 

J — OO 


2 


(B4) 
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C. Harmonics 


1. Scalar Harmonics: 


A set of dimensionless covariant harmonics Q can be defined on any LRS background as the eigenfunctions of the 
2-dimensional Laplace-Beltrami operator [ 8 , 9]: 


6 2 Q = —yQ, Q = Q = 0 (fc 2 > 0 ). 


The radial parameter r can be covariantly defined by 


f = r = -Ms-^0], 6 a r = 0. 


Any first order scalar may now be decomposed into this scalar harmonic basis 

^ = E4 fe) Q (fe) = M- 


(Cl) 


(C2) 


(C3) 


2. Vector Harmonics: 


Even parity vector harmonics can be systematically constructed from the scalar harmonics 


Qi k) = r5 a Q( k) =► Q a = Q a = 0 , 

(C4) 

5 2 Qa = ^ (1 - k 2 ) Q a , 

(C5) 

and odd parity vector harmonics are defined by 


Qi k) =re ab 6 b QW ^Q a = Q a = 0, 

(C6) 

S 2 Qa = 4r (1 - k 2 ) Q a . 

(C7) 

By construction, these harmonics are orthogonal parity inversions of each other, Q a 
Q a Qa = 0 . A crucial difference arises in that the even parity harmonics are not 
harmonics are: 

— e abQ Qa — e abQ and 

solenoidal but the odd parity 

k 2 

6 a Q a = -Q, 

r 

(C8) 

S a Qa = 0 , 

(C9) 

e ab 6 a Q b = 0, 

(CIO) 

e ab 6 a Q b = —Q. 

r 

(Cll) 

Any first order vector may be decomposed into this basis 


1>a = Y, 4 fc) Q ( a ] + = i’vQa + jvQa- 

(C12) 


k 


3. Tensor Harmonics: 

Likewise, we can extend these definitions to construct even and odd parity tensor harmonics 

Q ab — ^{a^b}Q = t > Q ab — Qab — 0? 

Q ab — ^c{a^ ^6} Q Qab — Qab — 0. 


(C13) 

(C14) 


23 


As with the vector harmonics, these tensor harmonics will be orthogonal parity inversions of one another: Q ab = 
— Cc{aQf,} C 44 Qab = e c { aQi ,} C and Q ab Qab = 0. Any first order tensor may be harmonically expanded in this basis 

ab = J2 Qab + T } Qab = Qab + Qab ■ (C15) 

k 


4. Relation to Spherical Harmonics 

The scalar covariant harmonics are defined on a sphere of radius r, whereas the conventional spherical harmonics 
are defined on the unit sphere. This allows us to construct a correspondence between the covariant and spherical 
harmonics. 

For scalar harmonics: 

Q = Y. (C16) 

For vector harmonics: 

Qa tT a , 

Qa = rY a . 

For tensor harmonics: 

Qab = r 2 Y ab , (C19) 

Qab = - \r 2 Y ab , (C20) 

where the last line depends on the exact definition of the axial tensorial harmonics (i.e. sometimes the prefactor of 
— 1/2 is omitted). 


(C17) 

(CIS) 



